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ABSTRACT 


’  This  is  a  review  of  equations  available  for  correcting  Pq  ,  ,  and  pH 

for  temperature  and  for  estimating  blood-gas  and  acid-base  variables  from 
measured  values: 

(  1.^  Oxygen  saturation  (Sq  )  from  Pq  ,  P^q  ,  pH,  and  body  temperature 
(T);  -  - 


2.  oxygen  concentration  from  Sq  ,  Pq  ,  and  hemoglobin  concentration 

r  1^  1 1^ 


3*^  Pq  at  half-saturation  (Pjq)  from  Pq  »  Sq  ,  pH,  P^q  ,  and  T; 

<  ,2  ,  '  r2  2  ,2 

<1.  "in  vivo  P^q"  from  Pq  »  Sq  »  and  T; 

.  /  V  <•  .2  '  i2 

5. ^-  base  excess  (BE)  from  pH,  P^q  ,  and  Hb; 

6. ''  "m  vivo  BE"  (BE^)  from  pH  and  P^q  ; 

r  1  (  A  2 

7.  "compensated  BE,"  from  P^^  ;  and 

7  2 

S.  bicarbonate  and  carbon  dioxide  concentrations  from  pH  and  P^q  • 

,  2 

Physiologic  considerations  are  emphasized,  with  mathematical  background  when 
it  contributes  to  understanding.  Algorithms  are  variously  compared  with  their 
graphic  counterparts,  with  the  data  from  which  they  were  derived,  and  with  each 

I 

other,  r, 


( 


INTRODUCTION 


With  increasing  availability  of  computers  and  programmable  calculators, 
many  computations  previously  performed  in  the  physiology  or  blood-gas  labora¬ 
tory  with  analog  devices,  such  as  slide  rules  (Sever inghaus,  1966),  nomograms 
(Siggaard  Andersen,  1963;  Kelman  and  Nunn,  1966;  Severinghaus,  1976a),  and  line 
charts  (Bradley  et  al.,  1956;  Severinghaus,  1958  and  1964;  Kelman  and  Nunn, 
1966),  can  now  be  carried  out  more  conveniently  using  digital  techniques.  To  do 
this,  one  needs  numeric  algorithms  to  represent  certain  physical  and  chemical 
relationships  in  the  blood: 

1.  It  is  sometimes  useful  to  estimate  the  percentage  saturation  of 

hemoglobin  with  oxygen  (S^  ),  using  measured  P_  and  a  standard 

curve  describing  oxygen-hemoglobin  equilibrium.  Occasionally  it  is 

also  useful  to  estimate  from  S-.  .  Not  only  are  accurate 

^2  ^2 

numeric  representations  of  the  oxygen  hemoglobin  equilibrium  curve 
(OH EC)  needed  for  these  transformations,  but  also  methods  to 
account  for  shifts  in  position  of  the  OHEC  caused  by  temperature, 
pH,  and  Pj,^  . 

2.  Blood  gases  and  pH  are  usually  measured  with  electrodes 

I 

thermostatically  controlled  at  37°C.  Results  should  be  corrected  to 
the  patient's  body  temperature  in  order  to  describe  the  true  physio¬ 
logic  condition. 

3.  In  assessing  acid-base  status  from  blood  pH  and  P^^  ,  most 
cliniciar\s  and  physiologists  find  it  useful  to  have  a  number  represent¬ 
ing  the  metabolic  component  of  acid-base  derangement.  Base  excess 
(BE)  or  plasma  bicarbonate  concentration,  (HCOj'j,  is  often  derived 
for  this  purpose. 
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This  paper  is  to  present  algorithms  for  digitally  carrying  out  the  above 
corrections  and  transformations.  Also  included  are  evaluations  of  how  faithfully 
the  algorithms  represent  their  graphic  counterparts  or  the  data  from  which  they 
were  derived.  All  of  the  equations  are  sufficiently  simple  that  calculations  can 
be  carried  out  on  small  "hand-held"  programmable  calculators,  such  as  the 
Hewlett-Packard  HP-67  or  the  Texas  Instruments  Tl-59.  Alternatively,  appro¬ 
priate  algorithms  can  be  consolidated  into  a  comprehensive  computor  program 
for  calculating  or  correcting  large  numbers  of  blood-gas  and  acid-base  variables 
in  a  busy  laboratory. 

METHODS  AND  RESULTS 

I.  OXYGEN  HEMOGLOBIN  EQUILIBRIUM  CURVE  (OHEC) 

A.  The  Standard  OHEC 

The  sigmoid  curve  describing  relationships  between  P^  and  S^  in  blood  is 

usually  mathematically  modeled  from  3S  data  pairs  reported  by  Severinghaus 

(1%6)  (Table  1).  These  matching  values  of  S^  (1-99.95  percent)  and  P_  (1.9- 

^2  ^^2 

700  torr)  represent  the  OHEC  for  normal  adult  human  blood  at  37°C,  pH  =  7.<»0, 
=  40  torr.  In  table  2  are  nine  algorithms  for  calculating  Sq  from  Pq 
(Kelman,  1966;  Thomas,  1972;  Aberman  et  ^.,  1973;  Roughton  and  Severinghaus, 

4 

1973;  Lutz  et  ^.,  1975;  Ruiz,  1975;  Severinghaus,  1976b  and  1979),  and  in  table  3 

three  algorithms  for  calculating  P^  from  S^  (Lutz  et  al.,  1975;  Tien  and 

^2  ^2 

Gabel,  1977;  Severinghaus,  1979).  Equations  characterize  the  entire  OHEC, 
except  where  a  limited  range  is  indicated. 

Fidelity  with  which  an  algorithm  represents  the  data  from  which  it  was 
derived  is  estimated  from  residual  errors;  that  is,  differences  between  the 
original  data  and  values  calculated  from  the  equation.  In  table  4  is  the  range  of 
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residual  errors  resulting  when  each  equation  is  compared  with  the  Severinghaus 
OHEC.  Also  listed  for  each  equation  is  the  root  mean  square  of  residual  errors 
(RMS);  that  is,  the  square  root  of  the  mean  of  the  squares  of  all  residual  errors,  a 
measure  of  the  "average  error"  from  applying  the  algorithm.  Table  4  also 
contains  the  approximate  time  it  takes  to  compute  or  P-  with  each 

O2  Oj 

equation,  using  a  typical  small  programmable  calculator. 

Ruiz's  algorithm  for  estimating  from  ,  although  clearly  superior  to 
the  others  in  accuracy  (RMS  =  0.07  percent  saturation)  may  not  be  ideal  for  all 
applications  because  of  its  relatively  long  computation  time  (20  sec)  when  a 
small  calculator  is  used.  Equation  fil  of  Severinghaus  retains  reasonable 
accuracy  (RMS  =  0.34  percent)  while  offering  rapid  solution  (1.5  sec).  A 
compromise  represented  by  Aberman's  algorithm  (RMS  =  0.17  percent;  13.5  sec) 
may  be  preferred  by  some. 

For  estimating  P^.  from  S^.  ,  application  of  two  of  the  algorithms  in  table 

3,  Lutz's  low-range  equation  when  Sq  is  below  20  percent  and  Tien's  wide-range 

algorithm  when  Sq  is  20-97.5  percent,  gives  a  combined  root  mean  square  of 
2 

0.10  torr  and  residual  errors  from  -0.21  to  0.18  torr.  Although  residual  errors  for 

both  equations  are  low  where  they  overlap  between  20  and  40  percent  saturation, 

Tien's  algorithm  is  more  accurate  than  Lutz's  in  this  range. 

B.  Effects  of  Temperature,  pH,  and  P^q  on  Position  of  the  OHEC 

When  developing  algorithms  for  dealing  with  changes  in  the  affinity 

between  oxygen  and  hemoglobin  owing  to  temperature,  pH,  and  ,  it  is 

generally  assumed  that  these  variables  influence,  position  but  not  shape  of  the 

OHEC  (Severin^us,  1958,  1%4,  and  1966;  Roughton,  1964;  Astrup  et  a[.,  1965; 

Kelman,  1966;  Aberman  et  al.,  1975).  For  Sn  =  ^"95  percent  where  this 

assumption  is  valid,  P^  is  increased  or  decreased  by  a  single  multiplicative 

^2 
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factor  for  any  combination  of  temperature,  pH,  and  Pqq  deviating  from  the 

standard  37°C,  7.40,  and  40  torr.  Hence,  the  S^,  of  blood  having  nonstandard 

temperature,  pH,  and  can  be  approximated  from  the  standeu-d  OH  EC  by 

first  calculating,  using  an  equation  to  be  derived  in  the  next  three  paragraphs,  a 

virtual  Pq  (P^)  that  would  exist  if  temperature,  pH,  and  P^q  were  made 

standard  (Severinghaus,  1958  and  1964;  Kelman,  1966;  Kelman  and  Nunn,  1966; 

Thomas,  1972;  Ruiz,  1975).  The  value  of  P^  can  then  be  substituted  for  P  in  one 

of  the  algorithms  in  table  2  to  compute  the  expected  at  the  observed  values 

^2 

of  Pq  ,  P^Q  ,  pH,  and  temperature. 

When  Sp.  is  held  constant  and  the  affinity  between  oxygen  and  hemoglobin 
\  ^2 

is  modified  by  an  alteration  in  temperature  or  pH,  ratios  of  respective  values  for 
Pq  vary  as  the  power  (base  10)  of  the  change  in  temperature  or  pH  (Bradley  et 
al.,  1956;  Astrup  et  ^.,  1965;  Severinghaus,  1966;  Kelman  and  Nunn,  1966): 

P^  /P^^  =  expjQlffTj-Tj))  (1) 

PpH/PpH  =®’‘Pl0^8(pHi-pH2))  (2) 


where  P_  and  P,.  are  P^  "s  at  different  temperatures  (°C),  P  „  and  P  „  are 

*1  '2  ^2  P"l  P”2 

Pq  ’s  at  different  pH's,  and  f  emd  g  are  constants  derived  empirically.  Factors 
'^2 

reported  by  Severinghaus  (1966)  are  usually  used  for  these  conversions: 


f  =  AlogjpPQ  /AT  =0.024  (3) 

g  =  AlogjQPQ  /ApH  =  -0.48+(0.0013(BE)/ApH)  (4) 

Equation  4  acknowledges  that  changing  pH  by  altering  blood  [HCO^']  has 

less  influence  on  the  affinity  between  oxygen  and  hemoglobin  than  changing  pH 

by  altering  ^C02  (Naeraa  et  ^.,  1963  and  1966;,  Severinghaus,  1966).  Kelman 

(1%6),  in  his  equation  for  calculating  virtual  P^  ,  used  the  form: 

^2 

g*  =  AlogjjjPQ  /ApH  =  -0.40-(0.06(logjQ40-logjoPcO,^^^P”^ 
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To  produce  a  single  multiplicative  factor  accounting  for  effects  of  both 

temperature  and  pH  on  oxygen-hemoglobin  affinity,  equation  1  is  multiplied  by 

equation  2.  To  do  this,  the  exponents  f(Tj-T2)  and  g(pHj-pH2)  or  g'(pHj-pH2) 

are  added.  Then  standard  values  for  temperature  and  pH  are  introduced,  and 

values  for  f  and  g  or  g’  from  equations  3  and  4  or  4'  are  substituted.  This  gives 

two  algorithms  for  calculating  virtual  P-.  (P^  or  P'  )  from  measured  P^  (P  ); 

^2  ’  ’  ^2  ^ 


Pv  =  Pj^(expjQ(-0.024(T-37)+0.48(pH-7.4)-0.0013(BE))) 
Pv  =  P^(expjQ(-0.024(T-37)+0.40(pH-7.4)+ 
0.06(log,Q40-logjoP(3Q^))) 


(5) 


(50 


where  P  and  P  refer  to  virtual  P^  "s  calculated  from  combinations  of 
v  V  ^2 

equations  1,  2,  and  3  with  equations  4  and  4',  respectively. 

Values  for  the  dependent  variables  P,/P  and  P  /P  fall  within  one 

V  m  v  m 

percent  of  each  other  for  all  independent  variables  in  the  ranges:  pH  =  7.2-7.6, 

Pco  =  30-80  torr,  BE  =  -15  to  +10  mEq/L,  and  Hb  =  10-20  gm/100  ml.  The 

ratios  fall  within  1.5  percent  of  each  other  when  pH  =  7. 2-7. 6,  P^q  =  20-80 

torr,  BE  =  -15  to  +15  mEq/L,  and  Hb  =  10-20  gm/100  ml. 

II.  ANAEROBIC  TEMPERATURE  CORRECTIONS 

With  changes  in  temperature  come  alterations  in  both  solubility  of  gases 

and  dissociation  const£mts  of  ionizable  solvents  and  solutes.  Therefore,  Pp.  , 

^2 

^COj  ,  and  pH  of  blood  change  as  temperature  is  anaerobically  altered.  With 
cooling,  pH  is  increased,  P^q  and  Pq  decreased;  with  warming,  the  opposite  is 
true. 


A.  eh 

Rosenthal  (1948)  reported  that  pH  of  whole  blood  decreases  0.0147  as 
temperature  is  elevated  1°C;  that  is,  ApH/AT  =  -0.0147.  Adamsons  and 
colleagues  (1964)  subsequently  showed  that  A  pH/ AT  is  not  constant  but  is 
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dependent  upon  pH  and  carbon  dioxide  concentration  (C^q  )  of  the  blood, 
reported  that: 

ApH/AT  =-0.0146-0.005(pH3g-7.4)+(5xl0"^)(C^Q  -20) 
where  pH^g  is  pH  at  38®C.  Severinghaus  (1966)  approximated  this  with 
ApH/AT  =  -0.0146-0.0065(pH3g-7.4)+(3xl0'^KBE) 


They 

(6) 

(7) 


and  pointed  out  that  the  final  term  may  be  ignored,  because  the  error  in  pH 
correction  owing  to  abnormal  BE  would  be  only  +0.006  in  the  extreme  case  when 
AT  is  +10°C  and  BE  is  +20  mEq/L. 

The  correction  per  degree  centigrade  will  be  insignificantly  different  if 
temperature  of  the  pH  electrode  is  37  rather  than  38°C.  Therefore,  pH  at  T°C 
can  be  calculated  from: 

pH^  =  pH37-(T-37)(0.0I46+0.O065(pH37-7.4))  (8) 


When  blood  is  anaerobically  warmed  or  cooled,  remains  constant 

(Bradley  et  1956;  Severinghaus  et  ^.,  1956;  Severinghaus,  1966),  leaving 
changes  in  pH,  dissociation  constants,  and  solubility  of  carbon  dioxide  to  account 
for  temperature-dependent  alterations  in  .  As  with  ,  temperature 
correction  for  P^q  follows  an  exponential  relationship  (Bradley  et  al^.,  1956; 
Nunn  et  al.,  1965;  Kelman  and  Nunn,  1966;  Severinghaus,  1966): 


P'^  /P'j  =  exp^Q(f'(T^-T2))  (9) 

where  P_.  and  P_  are  P^,..  "s  at  different  temperatures. 

*1  *2  ^°2 

Although  Bradley,  Stupfel,  and  Severinghaus  (1956)  reported  theoretical 
values  for  f  ranging  from  0.0175  to  0.0222,  depending  upon  pH  and  temperature, 
Nunn  and  colleagues  (1965)  experimentally  found  that  f  =  0.019  for  wide  ranges 
of  Pco^  temperature  (18.0-36.5°C).  Severinghaus  (1966) 

subsequently  advocated  use  of  0.019  for  temperature  correction  near  37°C,  with 
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the  admonition  that  correction  factors  for  lower  body  temperatures  (less  than 
30°C)  will  be  too  large,  especially  if  blood  pH  is  cilso  low.  For  body 
temperatures  of  30-40°C,  the  following  algorithm  can  be  used  to  correct 
for  differences  in  temperature  between  the  measuring  electrode  and  the  body; 

fXj  =  F  ,^expjQ(0.019(T-37))  (10) 

where  P'y  is  P^q  at  body  temperature  (T°C)  and  P'^^  is  P^q  measured  at 
37°C. 

Bradley  and  colleagues  (1956)  calculated  values  for  P^/P'^^  using  the 
following  transformation  of  the  Henderson-Hasselbalch  equation: 

Pco^  =  Cco^Aa-fexpiolpH-pKl+l))  (11) 

where  a'  is  the  solubility  of  carbon  dioxide  in  plasma.  Because  remains 

constant  during  change  in  temperature  (Severinghaus  et  1956)  and  will  cancel 
out  in  numerator  and  denominator,  the  ratio  of  values  for  Pqq  at  T°C  and 
37°C,  expressed  as  in  equation  11,  is: 

=«'37(expjQ(pH37-pK37)+l)/  (12) 

a'3.(expjQ(pHy-pKj)-t-l) 

With  substitution  of  currently-accepted  values  for  a  and  PK37  (37°C)  and 
and  pKy  (T°C)  at  various  values  of  pH  (7. 0-7.6)  and  temperature  (35-40°C) 
(Severinghaus,  1964),  ratios  P'y/P'3^  arising  from  this  equation  fall  within  1 
percent  of  those  calculated  from  equation  10.  For  correcting  to  progressively 
lower  temperatures,  ratios  P'.j./P'3^  calculated  from  equation  10  become  sub¬ 
stantially  larger  than  those  predicted  from  the  Henderson-Hasselbalch  equation 
(equation  12),  especially  if  pH  is  low.  For  example,  when  correcting  P^q  of 
blood  h  ng  pH  =  7.0  from  37  to  30°C,  measured  P^q  would  be  multiplied  by 
0.736  acv.ording  to  the  algorithm  (equation  10),  by  0.712  according  to  the 
Henderson-Hasselbalch  equation,  a  difference  of  3.4  percent. 
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Change  in  the  Pq  ot  blood  with  anaerobic  cooling  or  warming  in  the 
2 

midrange  of  the  OH  EC  is  caused  primarily  by  an  alteration  in  the  affinity 
between  oxygen  and  hemoglobin;  that  is,  in  position  of  the  OHEC  (Severinghaus, 
1966).  For  Sq  up  to  about  75  percent,  the  Pq  correction  for  temperature  is 
dependent  on  changes  in  position  of  the  OHEC  with  temperature  (equation  3)  and 
pH  (equation  4)  and  on  changes  in  pH  with  temperature  (equation  7).  Multiplying 
equations  4  and  7  gives  AlogjgP^  /AT  ov.ing  to  change  in  pH  with  temperature: 

AlogjoPo^/AT  =0.007^.0031(7.4-pH37)-(lxl0'^XBE)  (13) 

Adding  this  to  equation  3  produces  an  algorithm  accounting  for  alterations  in 
Pq  caused  by  a  combination  of  the  effects  of  temperature-dependent  changes 
in  pH  and  direct  effects  of  temperature  on  hemoglobin-oxygen  affinity 
(Severinghaus,  1966): 

AlogjgPQ^/AT  =  0.031+0.0031(7 .4-^>H37)-(IxlO■^XBE)  (14) 

The  second  and  third  terms  of  this  equation,  having  to  do  with  dependence  of 
A  pH/  AT  on  pH  and  BE,  may  be  ignored,  since  they  would  produce  a  maximum 
change  in  AlogjpPQ  /  AT  of  only  _^.0018  (^.4  percent  error  in  the  ratio  of 
Pq^'*  different  temperatures,  P.j./P3^,  per  degree  centigrade)  for  variations  in 
pH  of  ^.5  and  in  BE  of  +20  mEq/L.  Hence,  for  midrange  of  the  OHEC 
(Severinghaus,  1966): 

A  log]  oPq^/ AT  =0.031  (15) 

This  assumes  constancy  not  only  of  oxygen  concentration  (C^  )  but  also  of 

Sq  during  anaerobic  temperature  change.  Saturation  greater  than  75-80 
2 

percent  (P^.  above  40-45  torr)  cannot,  however,  be  assumed  to  remain  constant 
^2 

during  change  in  temperature.  When  blood  is  cooled,  some  dissolved  oxygen 
combines  with  hemoglobin,  because  during  cooling  the  increase  in  affinity 
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between  hemoglobin  and  oxygen  is  dominant  over  the  increase  in  solubility  of 

oxygen  in  blood  (Hedley-Whyte  and  Laver,  1964;  Nunn  et  1965;  Severinghaus, 

1966;  Roughton  and  Severinghaus,  1973).  The  resulting  slight  elevation  in  S~ 

becomes  more  influential  near  the  top  of  the  OHEC,  where  small  increases  in 

Sq^  result  in  large  increases  in  Pq  .  Therefore,  at  Pq  values  greater  than  40- 

45  torr,  a  reduction  in  temperature  produces  proportionally  less  decrease  in  P^ 

fig.  la,b  than  at  iower  Pq  's  (fig.  la,b).  That  is,  the  reduction  in  P^  with  cooling,  caused 

2  ^^2 
by  both  increased  solubility  of  oxygen  and  increas^  oxygen-hemoglobin  affinity, 

is  attenuated  at  higher  Pq  's  by  the  augmentation  of  induced  by  cooling. 

2  ^^2 

The  same  is  true  for  warming,  during  which  the  increase  in  P_  is  attenuated  at 

2 

higher  Pq  's  owing  to  the  transfer  of  small  amounts  of  oxygen  from  the 
^2 

combined  to  the  dissolved  state. 

Variable  increasing  Pq  and  Sq  is  accounted  for  in 

the  data  of  Severinghaus  (1966)  that  underlie  the  edgorithm  reported  by  Ruiz  and 

colleagues  (1975)  for  calculating  P_  at  T°C  (Pj  from  P^  measured  at  37°C 

2  ^2 

(P37): 

P^  =  P37(expjQ((T-37)(0.0265/((P37/146)^+l)+  (16) 

0.0007/(0.02(P37-230)^+l)+0.0047))) 

-  An  edgorithm  using  S^.  (S)  rather  than  P^  comes  from  a  commercial  source 

^2  °2 

(HP-67/HP-97  Users'  Library  Solutions.  Pulmonary.  Hewlett-Packard,  1000  N.E. 
Circle  Boulevard,  CorvaJis,  OR  97330): 

P^  =  p37(expjQ((T-37X3130-62.5(S)+0.312008(S^))/  (17) 

(100,000-1,993(S)+9.9313{S^)))) 

The  terms  accounting  for  variable  blog^pP^  /LI  in  equations  16  and  17 

agree  with  each  other  to  within  2  percent  when  any  of  the  first  35  pairs  of  values 

for  Pp.  (1.9-225  .torr)  and  Sp.  (1-99.5  percent)  on  the  standard  OHEC  (table  1)  is 
U2  O2 
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substituted  in  them.  More  useful  for  comparison  than  percentage  differences  in 
^IoSjqPq  /at,  however,  are  percentage  differences  in  which  are 

dependent  upon  T-37  and  which  can  be  calculated  after  exponential  transfor¬ 
mation  of  (T-37KAIogjQPQ  /AT); 

P^  =P37(expjQ((T-37)(AJogjQPQ  /AT)))  (18a) 

=expjQ((T-37XAlogjQPQ^/AT))  (18b) 

When  Sq  =  1-99.5  percent,  factors  PJ/P37  computed  from  equations  16  and  17 
differ  by  less  than  0.1  percent  per  degree  centigrade  change  in  temperature. 

For  Sq  up  to  99  percent  (Pq  to  159  torr),  values  of  Alog^pPQ  /  AT 
calculated  from  equations  16  and  17  agree  to  within  2  percent  with  those  read 
from  the  Severinghaus  graph  of  AlogjQPQ  /  AT  plotted  against  Sq  and  Pq 
(Severinghaus,  1966)  (fig.  la).  As  above,  it  follows  that  factors  P-J-/P37  from  the 
equations  differ  with  those  from  the  graph  by  less  than  0.1  percent  per  degree 
centigrade  change  in.  temperature.  Therefore,  either  equation  16  or  17  Cein  be 
used  to  represent  the  factors  recommended  by  Severinghaus  in  1966  to  correct 
P-.  for  anaerobic  temperature  change. 

Severinghaus  (1979)  has  recently  published  his  own  algorithm  for  this 
correction: 

P.J.  =  P37(expg((T-37K0.013+0.058/(U0.243(P37/100)^-**))))  (19) 

(The  exponent  here  is  to  the  base  e,  a  fact  that  will  be  accounted  for  later,  when 

this  equation  is  compared  with  algorithms  having  exponents  to  the  base  10.)  This 

equation  was  derived  theoretically  from  an  expression  for  slope  of  the  OH  EC  as 

a  function  of  Pq  .  Correction  factors  P‘j-/P37  calculated  from  equation  19  are 

lower  by  as  much  as  0.36  percent  per  degree  centigrade  (at  P-.  =  225  torr,  S^ 

^2  ^^2 

=  99.5  percent)  and  higher  by  up  to  0.28  percent  per  degree  centigrade  (P^'  = 

^2 

99.6  torr,  S^  =  97.5  percent)  than  those  calculated  from  equation  16  or  17.  The 
'^2 
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fls.  Ic 


three  algorithms  might  be  considered  comparable  when  accounting  for  small 
temperature  variations. 

From  data  of  Nunn  and  colleagues  (1965),  Kelman  and  Nunn  (1966) 
developed  another  algorithm  for  calculating  at  a  temperature  other  than 
that  at  which  P^.  is  measured: 

P^  =  P37(expjQ((T-37)(0.0052+0.0268(l-exp^(-0.3(100-S))))))  (20) 

Values  for  Alog^gPQ  /  AT  =  0.0052+0.0268(l-exp^(-0.3(100-S)))  are  as  much  as  27 

percent  lower  than  those  resulting  from  applying  the  appropriate  segments  of 

equations  16  and  17  or  from  reading  values  from  the  Severinghaus  graph  (1966) 

(fig.  la);  values  of  Alog^gP^  /AT  from  equation  20  are  up  to  29  percent  lower 

than  those  from  equation  19.  However,  the  exponential  transformation  results  in 

values  of  Pj/P^^  from  equation  20  that  differ  by  no  more  than  1.3  percent  per 

degree  centigrade  (T-37)  with  those  from  equation  16  or  17,  and  by  a  maximum 

of  1.6  percent  per  degree  centigrade  with  those  from  equation  19.  Because 

percentage  differences  between  the  Kelman-Nunn  and  the  Severinghaus  values 

for  AlogjQPQ  /AT  are  maximum  when  Pq  =  110-130  torr  (fig.  la),  percentage 

differences  between  the  respective  correction  factors  increase  to  more 

than  5  percent  for  temperature  variation  (T-37)  greater  than  4°C  in  this  range  of 

P^  (fig.  Ic).  Considerations  for  choosing  between  algorithms  derived  from  the 
2 

Severinghaus  (1966)  or  the  Nunn  data  (1965)  are  beyond  the  scope  of  this  paper. 
III.  BASE  EXCESS  (BE) 


Several  algorithms  for  calculating  BE  with  computers  are  available  (Dell  et 
al.,  1967;  Maas  et  1972;  Thomas,  1972).  Two  others  can  be  solved  using  small 
programmable  calculators.  Siggaard  Andersen's  algorithm  (1966)  is  the  simpler 
of  the  two: 


BE  =  (l-0.0143(Hb))((HCO3“)-(9.5+1.63(Hb))(7.4-pH37)-24) 


(21) 


12 


fig.  2 


where  Hb  is  the  concentration  of  hemoglobin  (gm/lOOml).  Substituting  a  form  of 
the  Henderson-Hasselbalch  equation  to  express  [HCO^~]  in  terms  of  pH  (pH^^) 


and  ^***37^’  measured  at  37°C: 


BE  =(l-0.0143(Hb))(0.0306(P'j7KexpjQ(pH37-6.l))-  (22) 

(9.5+ 1 .6  3(Hb)K7 .4-pH37)-24) 

Included  are  constant  values  for  apparent  pK  (6.1)  and  solubility  of  carbon 
dioxide  in  plasma  (0.0306  mM/(liter  x  torr)).  No  improvement  in  accuracy  results 
from  correcting  pK'  lor  pH  using  a  linear  regression  of  data  reported  by 
Sever inghaus  (1964);  pK'37  =  6.4 104-0.042 14(pH37). 

Values  of  BE  calculated  with  Siggaard  Andersen's  algorithm  are  generally 
more  positive  or  less  negative  than  those  read  from  the  Severinghaus  Blood  Gas 
Calculator  (BGC)  (1966)  .  For  the  following  ranges:  BE  =  -15  to  +15  mEq/L, 
=  20-80  torr,  pH  r  7.2-7.6,  and  Hb  =  10-20  gm/lOOmI,  equation  22 
underestimates  negative  BE  by  up  to  1.9  mEq/L  and  overestimates  positive  BE 
by  as  much  as  1.3  mEq/L^.  Percentage  error  is  greatest  at  BE  =  0  mEq/L,  where 
overestimation  is  1.4  mEq/L  when  P^q  =  21.1  torr,  pH  =  7.60,  and  Hb  =  10 
gm/iOOml  (table  5). 

Gershwin  and  colleagues  (1974)  developed  an  iterative  digital  technique  for 
calculating  BE  from  the  point  of  intersection  between  a  mathematically-modeled 

I 

blood  buffer  line  and  the  base  excess  curve  on  Siggaard  Andersen's  acid-base 
nomogram  (1962)  (fig.  2).  Slope  of  the  buffer  line  is  approximated  by: 

s  =-1.16(exp^(0.0204(Hb)-0.01434(BE)))  (23a) 

pH  at  which  the  buffer  line  intersects  the  base  excess  curve  (fig.  2)  is  fomd  by 
simultaneously  solving  equations  representing  the  two  lines,  each  expressed  as  a 


function  of  Alog|QP^Q^/  and  pH; 
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pH'  =  (-B-(B^-4AC)°*^)/2A 


where 


A  =  s+0.77 


B  =logjQp^Q  -s(pH)-7.55-7.165(s4-0.77) 

C  =-7.165(logjQP^Q  )+7.165(s)(pH)+7.55(7.165)+0.058 
To  calculate  BE,  the  value  of  pH'  is  substituted  in  an  empiric  equation  for  the 
base  excess  curve,  expressed  as  a  function  of  pH  alone: 

BE  =  59.4(logjQ(3.36(pH'-0.1)))  (23c) 

In  practice,  an  initial  estimate  of  0  mEq/L  is  substituted  for  BE  in  equation 
23a.  Subsequent  estimates  obtained  from  23c  are  then  reintroduced  into  23a 
until  equations  23a,  b,  and  c  have  been  sequentially  applied  three  times.  (These 
calculations  take  about  18  sec  on  the  HP-67  programmable  calculator.)  The 
result  after  three  iterations  is  considered  the  best  estimate  of  BE;  results  from 
further  iterations  sometimes  diverge  from  the  real  value. 

Gershwin's  algorithm  probably  provides  the  best  available  method  for 
calculating  BE  using  small  programmable  calculators.  When  BE  =  -15  to  +10 
mEq/L,  P^^q  =  20-80  torr,  pH  =  7.2-7.6,  and  Hb  r  10-20  gm/lOOml,  it  gives 
values  for  BE  that  are  within  +0.8  mEq/L  of  values  read  from  the  BGC  (table  5). 
For  perspective,  differences  in  BE  of  +0.8  mEq/L  can  be  produced  by  simul- 
,  taneous  variations  in  P_^  and  pH  of  +0.5  torr  and  +0.005,  respectively,  at 

extremes  of  the  above  ranges. 

If  Hb  is  low  and  **C02  high,  equations  23a,  b,  and  c  consistently  overesti¬ 
mate  values  of  BE  above  +3  mEq/L.  Errors  become  progressively  greater  as  BE 
fig.  3  increases  and  can  become  considerable  when  BE  exceeds  +10  mEq/L  (fig.  3). 

When  Hb  =  10  gm/lOOmI  and  P^q  =  80  torr,  the  overestimation  is  2A  mEq/L  at 
BE  =  +15  mEq/L,  pH  =  7.347,  and  is  4.1  mEq/L  at  BE  =  +20  mEq/L,  pH  =  7.393. 
For  Hb  =  5  gm/lOOml,  errors  at  the  above  P^q  and  BE  are  +4.4  and  +8.8  mEq/L, 


14 


respectively.  When  Hb  =  5  gm/lOOml  and  s  60  torr,  the  algorithm 

overestimates  BE  by  2.6  mEq/L  at  BE  =  +15  mEq/L,  pH  =  7.433,  and  by  5.4 
mEq/L  at  BE  =  +20  mEq/L,  pH  r  7.4S0.  Therefore,  Gershwin's  algorithm, 
although  quite  good  otherwise,  should  not  be  used  for  severely  anemic  patients 
with  substantial  degrees  of  partially-compensated  respiratory  acidosis  or  meta¬ 
bolic  alkalosis. 


There  are  limitations  to  the  usefulness  of  BE,  however  determined,  for 
planning  treatment  of  patients  with  metabolic  acid-base  derangements.  A  major 
problem  is  that  BE  is  defined  according  to  the  way  pH  of  fully-oxygenated  blood 
changes  as  is  systematically  altered  in  vitro  (Astrup  et  ^.,  1960).  Many 

have  pointed  out  that  m  vivo,  where  carbon  dioxide  and  bicarbonate  are 
distributed  in  interstitial  fluid  as  well  as  in  blood,  BE  exhibits  a  different 
relationship  to  pH  and  P^q  than  in  the  test  tube  (Schwartz  and  Reiman,  1963; 
Brackett  et  al.,  1965  and  1969;  Brown  and  Clancy,  1965;  Michel  et  ^•,  1966;  van 
Ypersele  de  Strihon  et  al.,  1966;  Arbus  et  ^.,  1969).  Severinghaus  (1976a) 
proposed  that  the  "in  vivo  BE"  (BE^)  be  read  from  Siggaard  Andersen's  alignment 
nomogram  at  3  gm/lOOmI,  an  "average"  hemoglobin  concentration  in  the 
extracellular  fluid.  Hence,  when  bicarbonate  and  carbon  dioxide  re-equilibrate 


and  redistribute  as  P^q  ;  in  the  blood  is  altered,  the  BE^ 


remains  relatively 


constant,  signifying  constancy  of  the  metabolic  acid-base  state.  The  convention¬ 


al  BE  would  change  under  these  circumstances,  deceptively  implying  an  increase 
or  decrease  in  nonvolatile  acid. 


Severinghaus'  algorithm  (1976a)  for  approximating  BE^  is: 

BE3  =  37(expg((pH-7.4+0.345(log^(P^Q  /40)))/ 
(0.55-0.09(logg(Pco2^'^°^^^^' 
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For  BE  =  -20  to  +15  mEq/L  and  P^q  =  10-80  torr,  BE^  calculated  with  this 
equation  is  within  +0.3  mEq/L  of  that  read  from  a  special  modification  of 
Siggaard  Andersen's  alignment  nomogram  (Severinghaus,  1976a)  (table  5). 
Differences  between  the  eUgorithm  and  the  nomogram  are  less  than  1  mEq/L  for 
BE  =  -30  to  +25  mEq/L,  P^q  =  10-80  torr. 

Severinghaus  (1976a)  also  reported  an  equation  for  calculating  the 
"compensated  BE^"  (BE^)  that  would  be  expected,  according  to  the  data  of 
Brackett  and  colleagues  (1%9),  after  full  renal  compensation  for  prolonged 
hypercapnia: 


BE^  =  15(log^(P^Q  /40))+0.095(P^q^-40)  (25) 

VThis  equation,  too,  matches  the  nomogram  to  within  +0.3  mEq/L  for  P^q  up  to 
80  torr.  Dividing  equation  24  by  equation  25  gives  an  indication  of  the  degree  of 
metabolic  compensation  for  respiratory  acidosis. 

For  those  who  find  the  plasma  bicarbonate  concentration  helpful  when 
clinically  managing  acid-base  disorders,  the  following  version  of  the  Henderson- 
Hasselbalch  equation  can  be  used: 


(  HC03"  ]  =a'(P^Q^XexpjQ(pH-pK'))  (26) 

where  pK'  anda'  (solubility  of  carbon  dioxide  in  blood  plasma)  are  appropriate  for 
temperature  and  pH.  For  37°C,  the  following  includes  a  value  fora'  and  an 
estimate  of  pK'  corrected  for  pH,  according  to  data  reported  by  Severinghaus 
(1%4): 


[  HC03']  =  0.0306(P^q  XexpjQ(1.04214(pH)-6.4104))  (27) 

Likewise,  carbon  dioxide  concentration  can  be  calculated  from  the  follow¬ 
ing: 

^C02  <28) 

where,  as  above,  pK'  anda'  are  appropriate  for  temperature  and  pH.  For  37°C; 


^€©2  =0.0306(P^Q^Xl+exp.Q(1.04214(pH)-6.4104))  (29) 
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DISCUSSION 

The  accuracy  with  which  S^  can  be  estimated  from  P^.  using  an 
algorithm  depends  not  only  on  how  faithfully  the  algorithm  represents  the 
standard  OHEC,  but  also  on  how  "standard"  the  blood  in  question  is.  Shape  and 
position  of  the  OHEC  can  be  quite  nonstandard  when  there  is  carboxyhemogiobin, 
methemoglobin,  or  an  abnormal  hemoglobin  present  (Darling  and  Roughton,  1942; 
Roughton,  1964  and  1970;  Enoki^^.,  1969;  Stamatoyannopoulos  et  al.,  1971). 

If  5  percent  of  the  hemoglobin  is  liganded  with  carbon  monoxide,  as  is 
common  among  urban  cigarette  smokers  (Stewart  et  al.,  1973),  ^50  of  an 

otherwise  standard  OHEC  is  reduced  from  26.6  to  25.1  torr,  and  the  curve 

V 

becomes  more  hyp>erbolic  and  less  sigmoid.  For  such  blood,  S^.  estimated  from 

2 

PZ  according  to  the  standard  OHEC  can  be  as  much  as  5.0  percent  (saturation) 
^2 

lower  than  the  true  saturation  of  available  hemoglobin  (that  not  combined  with 

carbon  monoxide).  For  example,  if  P^  is  17.7  torr  and  carbon  monoxide 

^2 

saturates  5  percent  of  the  total  hemoglobin,  oxygen  will  actually  saturate  31.0 

percent  of  the  remaining  hemoglobin,  whereas  estimated  from  the  standard 

2 

OHEC  will  be  26.0  percent. 

Likewise,  the  presence  of  methemoglobin  increases  the  affinity  between 
oxygen  and  hemoglobin  and  produces  a  more  hyperbolic  OHEC  (Darling  and 
Roughton,  1942;  Enoki  et  al.,  1969).  P^q  of  blood  containing  abnormal 
hemoglobins  can  be  as  low  as  4.6  torr  (hemoglobin  Hiroshima)  or  as  high  as  70 
torr  (hemoglobin  Kansas),  and  shapes  of  the  OHEC's  for  abnormal  hemoglobins 
arc  often  quite  different  from  the  shape  of  the  standard  curve  (Stamato- 
yannopoulas  et  ^.,  1971).  Furthermore,  position  of  the  OHEC  is  strongly 
dependent  upon  the  intraerythrocyte  concentrations  of  enzymes  (hexokinase, 
pyruvate  kinase)  and  substrates  (2,3-diphosphoglycerate  (2,3-DPG),  adenosine 
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triphosphate),  which  can  vary  substantially  in  many  pathologic  and  physiologic 
conditions  (Thomas  et  1974).  For  all  these  reasons,  estimated  from 
using  the  standard  OH  EC  should  be  interpreted  with  caution. 

Some  would  estimate  even  beyond  to  oxygen  concentration  ). 

^2  O2 

can  be  calculated  from  a  combination  of  measured  P-  ,  measured 
^2  O2 

hemoglobin  concentration  (Hb),  estimated  (from  virtual  P^  ,  as  explained  in 

^2  ^2 

section  IB),  assumed  solubility  of  oxygen  in  blood  (a)  (Hedley- Whyte  and  Laver, 
1964),  and  assumed  oxygen-carrying  capacity  of  hemoglobin  (g)  (Dijkhuizen  et 
al.,  1977); 

Co  =a(Pn  )+3(So  KHb)  (30) 

V  °2  °2  °2 

(The  following  values  are  conventionally  used:  a  =  0.003  mlO2/(100  ml  blood  x 

torr)  and  S=  1.34-1.36  ml02/gmHb.)  Besides  all  the  potential  errors  associated 

with  estimating  Sp,  from  P^,  ,  uncertainty  regarding  the  functional  value  of  a 
2  ^2 

for  a  given  blood  sample  (Theye,  1970;  Dijkhuizen  et  ^.,  1977)  makes  the 

estimated  Cp.  tenuous,  at  best. 

^2 

Before  calculating  virtual  Pq  from  equation  5  or  5'  when  estimating  Sp, 

2  ^2 
from  Pq^,  values  for  pH,  P^q  ,  and  Pq  should  be  corrected  to  body 

temperature  with  equations  8,  10,  and  16,  17,  19,  or  20,  respectively.  Contrary 

to  the  recommendation  of  some  (Ruiz  ^  ^.,  1975),  it  is  incorrect  to  calculate 

virtual  Pq^  from  values  of  pH,  P^q  ,  and  Pq  as  measured  at  37°C,  using  an 

equation  that  does  not  adjust  for  abnormal  body  temperature.  This  approach 

falseiy  assumes  that  physicaiiy  warming  or  cooling  the  blood  from  body 

temperature  to  37°C  in  a  blood-gas  analyzer  results  in  the  same  change  in  Pp. 

^2 

that  would  be  expected  from  shift  in  position  of  the  OHEC.  As  explained  in 

section  IIC,  Sq  actuaily  increases  when  blood  is  cooled  and  decreases  when 
2 

blood  is  warmed,  owing  to  transfer  of  oxygen  between  hemoglobin  and  the 
dissolved  state.  Therefore,  when  temperature  is  adjusted  to  37°C  for 
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measurement,  is  incremented  or  decremented  in  such  a  way  that  the  value 
of  Pp.  at  37°C  cannot  be  directly  transformed,  through  a  representation  of  the 
OH  EC,  to  at  body  temperature.  Rather,  the  P^  at  body  temperature  must 

O2  O2 

first  be  estimated  using  an  algorithm,  such  as  equation  16,  17,  19,  or  20,  th?:t 

corrects  for  change  in  Sp.  as  well  as  for  shift  in  position  of  the  OH  EC  with 

^2 

warming  or  cooling.  The  value  of  P*  at  body  temperature  could  then  be 

^2 

substituted  in  an  equation  defining  the  OHEC  under  in  vivo  conditions  of 
temperature,  pH,  and  P^q  .  In  practice,  however,  it  is  more  convenient  to 
calculate,  with  equation  5  or  5',  a  virtual  P_  expected  under  standard 
conditions  and  then  to  compute  from  an  algorithm  representing  the  OhEC  at 
37%  pH  =  7.40,  and  P^.^  =  40  torr  (table  2). 

In  contrast,  values  for  pH  and  P^q  as  measured  at  37°C  can  be  directly 
substituted  in  equations  22-25  when  estimating  BE,  BE^,  and  BE^.  Base  excess  is 
defined  in  terms  of  relationships  among  pH,  P^q  ,  and  hemoglobin  concentration 
in  fully  oxygenated  blood  at  38°C  (Astrup  et  ^.,  1960;  Siggaard  Andersen,  1962 
and  1963).  However,  BE  remains  constant  during  anaerobic  temperature  change, 
because,  despite  thermally-induced  alterations  in  pH  and  P^q  ,  the  concentra¬ 
tion  of  nonvolatile  acid  does  not  change  with  warming  or  cooling  (Severinghaus 
et  al.,  1956).  Therefore,  when  estimating  BE,  it  should  be  permissable  to  use 
values  of  pH  and  P^^  at  37°C  or  at  body  temperature.  For  reasons  both 
practiced  and  theoretical,  it  is  preferable  to  use  values  as  measured  at  37%, 
near  the  temperature  at  which  titrations  were  carried  out  for  experimental 
definition  of  BE  in  terms  of  pH,  P^q  ,  and  hemoglobin  concentration  (Astrup  et 
al.,  1960;  Siggaard  Andersen,  1962  and  1963). 

Algorithms  for  calculating  Pq  from  Sq  can  be  used  to  estimate  the  P^q 
of  a  blood  sample  from  one  set  of  measurements  of  Pq  ,  pH,  Pqq  ,  and  Sq  .  To 
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do  this,  measured  Pq  is  standardized  to  37°C,  pH  =  7.40,  and  =  40  torr, 

using  an  equation  for  calculating  virtual  P^.,  : 

^2 


P  =  P  (exp,.(-0.024(T-37)+0.40(pH-7.4)+ 
s  m  ‘^10 

0.06(logjQ40-logjQP^Q^))) 


(31) 


where  P  is  the  standardized  and  P„  the  measured  P^  .  As  when  calculating  P* . 
s  m  O2  V 

with  equation  5',  independent  variables  should  be  corrected  to  body  temperature 

before  substitution  in  equation  31.  The  Pq  that  would  be  expected  (P^)  from 

the  measured  ,  if  the  blood  under  study  had  standard  oxygen-hemoglobin 

equilibrium,  is  calculated  using  one  of  the  algorithms  in  table  3.  The  ratio  of 

observed-to-expected  P^,  (P,/P  )  provides  a  factor  by  which  all  P^,  *s  in  the 
v  ^2  s  e  KJ2 

midportion  of  the  OHEC  are  shifted  to  the  right  or  left  of  the  standard  curve. 
Therefore,  multiplying  P^q  of  the  standard  OHEC  (26.6  torr)  by  this  ratio  gives 
an  estimate  of  the  P^q  of  the  blood  under  standard  conditions  of  temperature, 
pH,  and  P^q  (Weiskopf  and  Severinghaus,  1972;  Aberman  1975). 

Aberman  and  colleagues  found  good  reproducibility  when  using  this  tech¬ 
nique  to  determine  the  P^q  of  92  blood  samples  having  Sq  between  20  and  90 
percent;  for  3-6  replications  on  each  sample,  standard  deviation  of  P^q  was  +1.0 
torr.  They  also  found  good  agreement  between  the  "one-point"  technique  and  the 

more  traditional  method  of  plotting  three  or  more  values  of  against 

5  2 

corresponding  P^  ,  with  or  without  Hill's  transformation  (1910)  ,  and  reading 
2 

P^Q  from  the  curve  that  best  fits  the  points  (Aberman  et  al.,  1975). 

Of  more  physiologic  relevance  than  the  standardized  ^50  described  above 

is  the  PjQ  of  blood  under  in  vivo  conditions  of  temperature,  pH,  and  P^q  .  To 

estimate  the  in  vivo  P^q,  only  Pq  and  Sq  need  to  be  measured.  Because  the 

one-point  technique  for  estimating  P^q  is  more  accurate  when  Sq  is  between  20 

and  90  percent  (Aberman  et  ^.,  1975),  Sq  of  the  blood  sample  should  be 

^2 


tonometerd  to,  if  not  already  in,  that  range.  Measured  is  corrected  to  body 
temperature  using  equations  16,  17,  19,  or  20  to  give  the  observed  (P  ), 
while  expected  Pq  (P  )  is  calculated  from  measured  Sp.  with  an  algorithm 
from  table  3.  Then,  just  as  for  the  standardized  P^q,  m  vivo  P^q  is  calculated  by 
multiplying  26.6  torr  by  the  ratio  of  observed-to-expected  P^ 
reflecting  influences  of  body  temperature  and  blood  pH  and  P^q  on  the  affinity 
between  oxygen  and  hemoglobin,  the  in  vivo  P^q  provides  a  more  realistic 
appraisal  of  the  availability  of  oxygen  to  body  tissues  than  does  the  standardized 


It  is  becoming  increasingly  more  evident  that  relationships  among  P_  , 

^2 

Sq  ,  pH,  Pqq  I  Hb,  2,3-DPG,  and  temperature  are  more  complex  than 

previously  thought.  For  example,  not  only  is  AlogjgP^  /  AT  reduced 

progressively  as  is  increased  above  75-80  percent  (fig.  la),  but 

^2 

AlogjQPQ  /AT  is  reduced  further,  when  Sq  is  above  90  percent,  by  hypercapnia 
and  by  anemia  (Kelman,  1968).  In  blood  with  a  low  concentration  of  2,3-DPG, 
AlogjQPQ  /  ApH  is  lower  than  normal  when  pH  is  changed  with  nonvolatile  acid, 
but  higher  than  normal  when  pH  is  changed  with  carbon  dioxide  (Wranne  ^  al., 
1972).  Whether  the  concentration  of  2,3-DPG  is  low  or  normal,  Alog|QPQ  /ApH 
owing  to  change  in  Pqq  w  greater  at  lower  than  at  higher  Sq  (Hlastala  and 
Woodson,  1975).  When  BE  is  substamtially  reduced  (for  example,  to  -10  mEq/L), 
AlogjQPQ  /ApH  is  additionally  increased  at  all  saturations  (Adams  et  al.,1968). 
These  interactions,  2md  others  like  them,  may  require  consideration  in 
physiologic  investigations,  but  they  probably  don’t  have  to  be  accounted  for  in 
algorithms  used  for  clinical  purposes. 

Kelman's  (1966)  and  Thomas'  (1972)  algorithms  for  estimating  Sp^  from 

^2 

Pp.  involve  the  same  polynomial  equation: 

^^2 
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S  =  100(ajP+a2P^+a3P^>P^)/(a^>a5P>a^P^+a7P^+P^)  (32) 

Using  the  following  variation  when  programming  for  small  calculators  can  result 
in  calculation  times  that  are  one-third  to  one-half  as  long  as  when  the  above 
expression  is  straightforwardly  programmed: 

S  =  100((P(aj+P(a2+P(a3+P)))))/(a^+P(a3*P(ag+P(a2+P))))  (33) 

Whereas  it  takes  only  2  sec  to  solve  the  specially-:-..ctored  version  of  the 
Kelman- Thomas  equation  on  the  HP-67  programmable  calculator,  the  conven¬ 
tional  form  takes  about  4  sec  when  powers  are  generated  by  sequential 
multiplication  and  about  6  sec  when  the  y*  function  is  used.  Computation  times 
in  table  4  for  polynomial  equations  of  Thomas,  Adair,  Kelman,  and  Tien  reflect 
this  programming  strategy. 

While  calculating  the  P^.  at  body  temperature  from  that  measured  at 
37”C,  it  is  convenient  to  correct  for  a  common  measuring  artifact.  Because  of  a 
diffusion  gradient  for  oxygen  through  liquid  in  a  conventional  oxygen  analyzer, 

the  measured  P^  is  lower  than  the  true  blood  P_  (Polgar  and  Furster,  1960; 

^2  ^2 

Moran  et  ^.,  1966;  Rhodes  and  Moser,  1966;  Adams  et  ^.,  1968;  Severinghaus, 
1968;  Hulands  et  ^.,  1970;  Severinghaus  and  Bradley,  1971).  This  difference 
would  not  be  seen  if  the  blood  were  continually  stirred  during  analysis.  Blood 


is  usually  underestimated  by  about  4  percent  when  currently-available 


commercial  electrodes  are  calibrated  directly  with  gas.  This  error  can  be 


corrected  with  a  factor  determined  by  measuring  the  P^  in  air  and  in  water 

^^2 


tonometered  with  air  at  37°C. 


To  account  for  slower  diffusion  of  oxygen 


through  blood  tham  through  water,  the  difference  between  P^  in  air  and  in 

^2 

water  is  multiplied  by  1.8.  The  difference  is  then  divided  by  the  P_  in  aur  to 

^2 

give  a  factor  with  which  to  correct  blood  Pq^  readings  for  "stirring  artifact": 
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CAPTIONS  FOR  FIGURES 


Fig.  1.  Factors  to  correct  Pq  for  anaerobic  temperature  change  are 

2 

dependent  upon  P^,  . 


Values  of  Alog^QP^  /  AT  from  Kelman  and  Nunn’s  algorithm  are  as 
much  as  27  percent  lower  than  those  reported  by  Severinghaus  in 
1966.  When  Sq  is  99  percent  or  less,  values  of  Alog^gP^  /AT 
calculated  from  Ruiz's  or  Hewlett-Packard's  algorithm  fall  within  2 
percent  of  the  Severinghaus  graph. 


To  correct  P^  from  37*^C  to  body  temperature,  the  measured  value 
2 

of  P^  is  multiplied  by  the  calculated  ratio  of  P^  at  body 
*^2  ^2 
temperature  (P^)  to  that  at  37°C  (P^^).  PjlPyj  calculated  from 

Ruiz's  or  Hewlett-Packard's  algorithm  differ,  for  any  P.  ,  with 

°2 

values  calculated  from  the  Severinghaus  graph  by  less  than  0.4 
percent  per  degree  centigrade. 


Values  for  P^/P^^  calculated  from  the  Kelman-Nunn  algorithm  differ 
by  as  much  as  1.3  percent  per  degree  centigrade  with  those 
calculated  from  the  graph  of  Severinghaus  (or  from  the  algorithms  of 
Ruiz  and  Hewlett-Packard).  Percentage  differences  are  maximum  at 
Pq^  =  110-130  torr.  Kelman  and  Nunn's  corrections  for  temperature 
are  lower  than  those  of  Severinghaus  when  Pq  is  greater  than  about 
65  torr  and  are  higher  than  those  of  Severinghaus  elsewhere  (fig.  la). 
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Fig.  2. 


V 

^  / 


Fig.  3. 


When  correcting  a  value  of  greater  than  65  torr  for  hypothermia, 
the  ratio  calculated  from  the  Kelman-Nunn  algorithm  is 

larger  than  that  calculated  from  the  Severinghaus  data;  that  is, 
PJ/P37  is  closer  to  unity,  which  represents  a  lesser  correction  (fig. 
lb). 


To  calculate  BE  using  Gershwin's  algorithm,  slope  of  the  blood  buffer 
line  is  approximated  by  substituting  a  value  for  Hb  and  an  estimate 
for  BE  in  an  empirically-derived  equation  (23a).  Equations  represent¬ 
ing  the  blood  buffer  line  (rectilinear)  and  the  base  excess  curve 
(hyperbolic),  each  expressed  in  terms  of  pH  and  »  are 

simultaneously  solved  for  pH,  eliminating  the  terms.  The 

resulting  value,  pH',  is  the  abscissa  of  the  point  of  intersection 
between  the  blood  buffer  and  base  excess  lines.  This  value  is  then 
substituted  in  equation  23c  expressing  BE  in  terms  of  pH  alone. 
Equation  23c  mathematically  represents  reflection  of  the  base  excess 
curve  on  the  horizontal  axis.  Respective  values  for  BE  and  pH, 
according  to  this  equation,  can  be  read  from  points  at  which  vertical 
lines  intersect  the  b«ise  excess  curve  and  the  horizontal  axis. 


When  P^Q  is  high  and  hemoglobin  concentration  is  low,  Gershwin's 
algorithm  overestimates  values  of  BE  to  a  progressively  greater 
degree  as  BE  and  P^q  increase  and  as  Hb  decreases. 


TABLE  1 

Data  Defining  Severinghaus'  Standard  Oxygen  Hemoglobin 
Equilibrium  Curve  (OHEC) 
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Algorithms  for  Calculating  S-  (S)  from  Po.  (P)  According  to  the  Standard  Oxygen  Hemoglogbin  Equiiibrium  Curve 


TABLE  3 

Algorithms  for  Calculating  P-  (P)  from  Sp.  (S)  According  to  the  Standard  Oxygen  Hemoglobin  Equilibrium  Curve 
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Deviations  from  the  Standard  Curve  and  Computation  Times  for  Algorithms  Describing 
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Table  4  (Continued) 
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less  than  the  full  OHEC  are  in  parentheses. 


TABLE  5 


Ranges  of  Base  Excess  (BE  and  BEj)*  in  mEq/L,  Calculated  from 
Algorithms  within  the  Ranges:  =  20-80  Torr,  pH  =  7.2-7.6 


BE  or  BEj  from 

BE  from 

BE  from 

BEj  from 

Nomogram  or  Blood 

Siggaard  Andersen's 

Gershwin's 

Sever  inghaus' 

Gas  Calculator 

Algorithm* 

Algorithm* 

Algorithm  •• 

-15 

-13.1 

to 

-14. 

-14.6 

to 

-15.1 

-15.2 

to 

-15.2 

(♦I.*) 

to 
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to 
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-3.7 
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♦  10.0 

to 

♦  10.3 
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to 
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♦  15 
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to 

♦  16.3 

♦  14.8 

to 

♦  17.4 
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to 

♦  15.3 

(♦0.3 

to 

♦  1.3) 

(-0.2 

to 

♦2.4) 

(-0.2 

to 

♦0.3) 

•  Compared  with  the  Blood  Gas  Calculator  in  the  range  of  Hb  *  10-20  gm/ 100ml 
Compared  with  Siggaard  Andersen's  alignment  nomogram,  Hb  =  3  gm/ 100ml 
•••In  parentheses  are  the  ranges  of  differences  between  values  calculated 
from  the  algorithms  and  those  read  from  the  nomogram  or  the  Blood  Gas 
Calculator. 


FOOTNOTES 


*  Supported  in  part  by  Grant  GM  15904  from  the  National  Institutes  of  Health, 
United  States  Public  Health  Service 

2 

The  BGC,  a  sliderule  having  one  degree  of  freedom  more  than  a  nomogram,  was 
designed  by  Severinghaus,  in  collaboration  with  Siggaard  Andersen,  to  represent 
the  original  BE  titration  data  better  than  could  Siggaard  Andersen's  alignment 
Vnomogram.  (Personal  communication;  John  W.  Severinghaus,  M.D.)  Therefore, 
tlw  BGC  constitutes  a  higher  standard  for  comparison  than  the  alignment 
nomogram. 


^BE  =  -13.1  from  the  algorithm  compared  with  -15.0  mEq/L  from  the  BGC  when 
P^Q  =  20  torr,  pH  =  7.302,  and  Hb  =  20  gm/ 100ml. 

4 

BE  =  +16.3  from  the  algorithm  compared  with  +15.0  mEq/L  from  the  BGC  when 
P^2q  =  39.2  torr,  pH  =  7.60,  and  Hb  =  15  gm/lOOmJ. 

^Sigmoidicity  of  the  OHEC  is  transformed  to  near  rectilinearity  by  plotting 

“gainst  logioPo2* 


Soj  (percent) 
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